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ABSTRACT 

We report the results of a diffraction-limited, photometric variability study of the central 5" x 5" of the 
Galaxy conducted over the past 10 years using speckle imaging techniques on the W. M. Keck 1 10 m telescope. 
Within our limiting magnitude of mx < 16 mag for images made from a single night of data, we find a minimum 
of 15 K[2.2 /im]-band variable stars out of 131 monitored stars. While large populations of binaries have been 
posited to exist in this region, both to explain the presence of young stars in the vicinity of a black hole and 
because of the high stellar densities, only two binaries are identified in this study. First is the previously 
identified Ofpe/WN9 equal mass eclipsing binary star IRS 16SW, for which we measure an orbital period of 
19.448 ±0.002 days. In contrast to recent results, our data on IRS 16SW show an asymmetric phased light 
curve with a much steeper fall-time than rise-time, which may be due to tidal deformations caused by the 
proximity of the stars in their orbits. Second is the WC 9 Wolf-Rayet star IRS 29N; its observed photometric 
variation over a few year time-scale is likely due to episodic dust production in a binary system containing 
two windy stars. Our sample also includes 4 candidate Luminous Blue Variable (LBV) stars (IRS 16NE, 16C, 
16NW, 16SW). While 2 of them show variability, none show the characteristic of LBVs large increase or 
decrease in luminosity. However, our time baseline is too short to rule them out as LBVs. Nonetheless, the 
lack of evidence for these stars to be LBVs and their coexistence with a significant surrounding population of 
well established Wolf-Rayet stars is consistent with needing only a single recent starburst event at the Galactic 
center to account for all of the known, young, massive stars. Among the remaining variable stars, the majority 
are early-type stars and three are possibly variable due to line of sight extinction variations. For the 7 OB stars 
at the center of our field of view that have well-determined 3-dimensional orbits, we see no evidence of flares or 
dimming of their light, which limits the possibility of a cold, geometrically-thin inactive accretion disk around 
the supermassive black hole, Sgr A*. 

Subject headings: Galaxy: center — infrared: stars — stars: variables: other — 



1. INTRODUCTION 

The stellar cluster at the Galactic center (GC) presents 
a unique opportunity to study the evolution and properties 
of stars within the sphere of influe nce of a 3 — 4 x 1 6 Mp > 
supermassive black hole (SMBH) dGhez et alJl200l 120051: 
Schfad eTet alJ l2003). Photometric variability offers a useful 
approach to a number of outstanding questions regarding this 
stellar population which is compos ed of a mixture of old gi- 
ants and young, m a ssive s tars. dR rabbe et a . 1991, IT9951 



Blum et al.1 ll996ailR 120031: iFiger et all 120031: iPaumard et all 
2001L [2004b, 2006). For example, light curves can easily re- 



veal close binary stars, which are relevant in several ways to 
our understanding of stars at the Galactic center. 

First, binaries on radial orbits that are disrupted by the 
central black hole may provide a mechanism for captur- 
ing young stars from large galacto-centric radii, where the 
conditions are conducive to star formation, and retaining 
them at the smaller less hospitable radii where m any young 
stars are found today dGould & Ouillenl [2003). Second, 
binary companions may facilitate the production of dust 
around the WC sub-class of Wolf-Rayet stars, which are 
massive post-main sequence stars undergoing rapid mass 
loss. While conditions in the hostile environment (high 
temperatures in particular) of the stellar winds do not fa- 
vor the formation of dust (William s et al.l 1 1987b . com- 
pression within wind-colliding binary systems could over- 
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come this challenge JWhite & Becked 19951: 1 Veen et al.l l 1998b 
IWilliams & van der Huchlll2000t iLefevre et alJl2005lT Third. 
binaries provide a direct measurement of stellar masses. This 
is especially helpful for the most massive stars in the Galactic 
center as it would assist our understanding of the recent star 
formation history. 

Another way in which a photometric variability study con- 
strains the recent star formation history, as well as our un- 
derstanding of massive star evolution, is the possibility of 
identifying luminous blue variables (LBVs). There are cur- 
rently only 12 confirmed Galactic LBVs and 23 additional 
candidates, with 6 can didates in the Gal actic center IRS 16 
cluster of stars alone (IClark et al.l 120051) . The LBV phase 
plays an important, although poorly constrained, role in stel- 
lar evolution, because, during this phase, stars experience sig- 
nificant mass loss, with rates of ~ 10~ 2 M Q /yr during erup- 
tions a nd as high as high as 10~ 4 S M w /yr during quie scent 
phases ([Abbott & Contill987tlHumphreys & Davidsonlll994t 
lMassevll2003l) . From a star formation history stand point, the 
LBV phase is notable because it is the first of several post 
main sequence phases that only the most massive stars (M > 
60-85M©) may go through before becom ing supernovae. 
Stars stay in this phase for only ~ 10 4 years dStothers & Chinl 
1996) before entering t he Wolf-Rayet phase, wh ich typically 
lasts a few x 10 6 years (Meynet & Maeder 2005). Less mas- 
sive stars (M > 40M Q ) will skip the LBV phase and become 
Wolf-Rayet stars, but on time-scales longer than that of the 
more massive stars that experienced an LBV phase. There- 
fore, in principle, the numbers of LBVs and WR stars can 
constrain recent star formation histories (e.g. Pauma rdet al.l 
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120061: lFigei1l2004l) . In this context the candidate LBVs at the 
Galactic center are perplexing in the context of the > 25 Wolf- 
Rayet stars located in their immediate vicinity, since in a sin- 
gle starburst event one would not expect to see any WR stars 
if the most massive stars are just now evolving through the 
LBV phase. This is similar to the problem p osed by the p res- 
ence of two LBVs in the Quintuplet cluster Figerj (|2004|) . If 
confirmed, the LBV candidates would suggest that this region 
has undergone multiple recent star forming events or that our 
understanding of LBV evolution is incomplete. 

The photometry of stars in close proximity to the SMBH 
can also be used to constrain the properties of a possible 
cold, geometrically-thin inactive accretion disk around Sgr 
A* which could explain the present-day low luminosit y of Sgr 
A* dNavakshin & SunvaevlBOOl ICuadra et al.ll2003l) . In the 
presence of such a disk, we would expect to see nearby stars 
eclipsed or reddened when they pass behind the disk. 

Ve ry few photometric va riable studies of the Galactic center 
exist. iTamura et alJ (0,999) introduced the idea that stars close 
to the Galactic center are expected to have a higher fraction 
of ellipsoidal 2 and eclipsing variable binaries than the stars 
in the solar neighborhood, but found very few variable stars 
and no binary stars. Se eing-limited studies (Tamu raet alJ 
119991: iBlum et al.1 [l996a) are limited to the brightest stars, 
due to stellar confusion caused by the high stellar densities 
and proper motions close to the central b lack hole. With high 
angular resolution data, lOtt et al.1 (1 19991) have identified the 
only known eclips i ng binary system in this region (see also 
iDePoy et all l2004t iMartins et all |2006|) . Furthermore, they 
suggest that as much as half of their sample (K <13) may 
be variable. However, the variability fraction decreases at 
smaller galactocentric radii starting from ~ 5", suggesting 
that even at a high resolution of 0."13 their sensitivity to vari- 
ability is limited by stellar confusion. 

In this paper, we present the results of a stellar variability 
study of the central 5"x5" of our Galaxy, based on ten years 
of K[2.2 fim] diffraction-limited images from the W. M. Keck 
I Telescope (9 = 0"05). The observations are described in §2, 
and the data and methodology to determine variability in §3. 
We discuss the variable star population in §4, which includes 
identification of asymmetries in the light curve of the eclips- 
ing binary star IRS 16SW and the discovery of a likely wind 
colliding binary star in IRS 29N, and summarize our major 
findings in §5. 



2. OBSERVATIONS 

K-band (A = 2.2/im, AA=0.4/im) speckle imaging obser- 
vations of the Galaxy's central stellar cluster were obtained 
with the W. M. Keck I 1 m telescope using the fa cility near- 
infrared camera, NIRC (Mat thews & Soiferl["l994l) . Observa- 
tions taken from 1995 to 2004 have been described in d etail 
elsewhere dGhez et al.lll998l 12000. 20051 iLu et alll2005l) and 
new observations on 2005 April 24-25 were conducted in a 
similar manner, resulting in diffraction-limited images. Each 
night several thousand short-exposure frames were taken in 
sets of ~ 200, with NIRC in its fine plate scale mode, which 
has a scale of 20.40 ± 0.04 mas pixel " 1 and a correspondin g 
field of view (FOV) of 5. "22 x 5. "22 (iMatthews et al.lfT996h . 
TableQ]lists the date and number of frames obtained for each 
of the 50 nights of observations used in this study. 

2 Ellipsoidal variables are non-eclipsing binaries that are elongated by mu- 
tual tidal forces ISterken & Jaschek 1 99611 . 



3. DATA ANALYSIS AND RESULTS 
3.1. Image Processing 

The individual frames are processed in two steps to create 
a final average image for each night of observation. First, 
the standard image reduction steps of sky subtraction, flat- 
fielding, bad pixel correction, optical distortion correction 3 , 
and pixel magnification by a factor of two are carried out 
on each frame. Second, the frames from each night of 
obser vation are com bined using the method of "Shift-and- 
Add" (IChristoull 19911) with the frame se lection and weighting 
scheme prescribed bv lHornsteinl (12006!). In short, each frame 
is analyzed for Strehl quality using the peak pixel value of IRS 
16C, and low quality frames, which do not improve the cu- 
mulative signal to noise ratio (SNR) for the observations from 
each night, are rejected. This typically leaves ^1600 frames 
for each night or 37% of the original data set (see column 3 
of Table [T}. The remaining frames for a given night are com- 
bined with Shift-and-Add in an average that is weighted by 
each frame's peak pixel value for IRS 16C. The final images 
have typical Strehl ratios of ^0.07 (see column 8 of TableQ}. 
The dataset from each night is also divided into three equiva- 
lent quality (and randomized in time) subsets to make three 
independent weighted Shift-and-Add image subsets, which 
are used to determine measurement uncertainties and to reject 
spurious sources. 

Sources are identified in individual images and cross- 
identified between ima ges us ing the strate gy de veloped in 
iGhez et al .1 (fl998l 12000112005b and lLu et al.1 (120061) . which for 
this study entails four separate steps. In the first step of the 
source identification process, we generate a conservative ini- 
tial list of sources for each night of data to help minimize spu- 
rious source detection. This is done usi ng the point-spread- 
function (PSF) fitting routine StarFinder (Di olaiti et al.ll2000T) 
to identify sources in both the average images and the sub- 
set images. StarFinder identifies sources through cross- 
correlation of each image with its PSF model, which, for our 
implementation, is generated from the two bright stars IRS 
16C and IRS 16NW. The initial source list for each night of 
data is composed of only sources detected in the average im- 
ages with correlation values above 0.8 and in all three subset 
images with correlation values above 0.6. In the second step 
of the source identification process, the source lists from all 
nights are cross-identified to produ ce a master list of sources, 
using a process that is described in lGhez et al.l (1 19981) and that 
also solves for the sources' proper motions. To further ensure 
that no spurious sources have been detected we require that 
sources be detected in a minimum of 13 nights 4 . 

Figure Q] displays the 131 sources contained in our final 
master list. In the third step of the source identification pro- 
cess, we return to the original images to search more aggres- 
sively for the sources on the master list that were missed in 
some of the images. We explicitly feed the master list of 
sources into StarFinder and search for only these sources at 
their predicted positions with more lenient criteria, which re- 
quire correlation values above 0.4 for both average and subset 
images. In the fourth and final step, we impose a restriction 
on our source detections to ensure photometric reliability: we 

3 http://www.keck.hawaii.edu/inst/nirc/Distortion.html 

4 The threshold for the minimum number of nights was chosen by looking 
for a drop in the distribution of the number of nights that the sources were 
detected in the first pass at source identification. A minor drop is seen at 13 
nights. The final results are not very sensitive to this choice and we therefore 
have made a fairly conservative choice. 
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TABLE 1 
List of Observations 



Date 


Frames" 


Frames b 


Num. Stars c 


Num. Stars d 


SNR e 


Strehl 




(Obs.) 


(Used) 


(Initial) 


(Final) 






1995 Jun 10 


1200 


425 


54 


66 


10.3 


0.08 


1995 Jun 11 


2700 


1604 


95 


110 


15.0 


0.06 


1995 Jun 12 


2100 


1082 


107 


108 


12.0 


0.04 


1996 Jun 26 


4200 


585 


1 16 


1 19 


19.8 


0.04 


1996 Jun 27 


2300 


1260 


117 


121 


20.5 


0.04 


1997 May 14 


3600 


1851 


63 


82 


16.7 


0.06 


1998 Apr 02 


2660 


1649 


119 


121 


16.3 


0.05 


1998 May 14 


4560 


1748 


96 


1 14 


16.7 


0.04 


1998 May 15 


7030 


1953 


41 


55 


9.9 


0.06 


1998 Jul 04 


2280 


943 


108 


1 14 


17.8 


0.08 


1998 Aug 04 


6270 


1469 


87 


107 


15.5 


0.05 


1998 Aug 05 


5700 


1592 


94 


103 


17.3 


0.07 


1998 Oct 09 


2660 


1188 


83 


99 


17.1 


0.08 


1998 Oct 11 


570 


450 


79 


96 


15.0 


0.05 


1999 May 02 


7030 


1589 


115 


1 16 


19.9 


0.09 


1999 May 03 


2090 


1264 


103 


1 14 


16.0 


0.07 


1999 Jul 24 


5510 


2239 


113 


119 


18.5 


0.11 


1999 Jul 25 


950 


788 


102 


1 14 


15.0 


0.06 


2000 Apr 21 


3040 


947 


90 


1 12 


18.5 


0.04 


2000 May 19 


9880 


1970 


64 


91 


10.2 


0.09 


2000 May 20 


7600 


2146 


81 


97 


16.5 


0.10 


2000 Jul 19 


8740 


1939 


111 


1 16 


17.1 


0.07 


2000 Jul 20 


3420 


1454 


111 


119 


19.8 


0.09 


2000 Oct 18 


2280 


1807 


82 


107 


15.1 


0.05 


2001 May 08 


1520 


889 


111 


124 


20.1 


0.04 


2001 May 09 


6270 


1990 


76 


100 


13.1 


0.08 


2001 Jul 28 


4180 


1752 


99 


104 


20.0 


0.13 


2001 Jul 29 


6080 


1751 


105 


115 


19.0 


0.07 


2002 Apr 23 


7410 


1669 


117 


119 


17.2 


0.05 


2002 Apr 24 


7790 


1882 


104 


113 


19.6 


0.06 


2002 May 23 


1900 


1249 


66 


84 


18.3 


0.07 


2002 May 24 


2660 


1537 


85 


100 


17.7 


0.09 


2002 May 28 


2850 


1866 


59 


80 


13.2 


0.06 


2002 May 29 


3420 


1552 


86 


104 


16.8 


0.07 


2002 Jun 01 


5510 


1992 


43 


53 


7.7 


0.09 


2002 Jul 19 


4370 


1115 


41 


57 


7.0 


0.07 


2002 Jul 20 


3990 


1355 


51 


64 


11.0 


0.06 


2003 Apr 21 


5130 


1799 


93 


107 


19.1 


0.04 


2003 Apr 23 


5320 


1970 


69 


87 


13.2 


0.05 


2003 Jul 22 


5130 


1718 


70 


94 


16.6 


0.08 


2003 Sep 07 


4560 


1795 


108 


112 


16.3 


0.07 


2003 Sep 08 


4370 


1223 


97 


110 


12.3 


0.07 


2004 Apr 29 


6840 


1181 


53 


68 


14.5 


0.11 


2004 Apr 30 


4180 


1203 


98 


105 


16.6 


0.05 


2004 Jul 25 


5320 


2007 


98 


110 


18.1 


0.08 


2004 Jul 26 


8550 


2309 


33 


38 


6.7 


0.08 


2004 Aug 29 


3230 


1328 


120 


122 


21.0 


0.10 


2005 Apr 24 


7410 


2195 


51 


60 


11.6 


0.07 


2005 Apr 25 


9500 


2035 


116 


116 


20.8 


0.05 


2005 Jul 26 


6650 


1497 


98 


113 


19.0 


0.06 



NOTE. — All observations are speckle K-band (A„ = 2.2/^m, AA=0.4^m) images. 
a The number of frames observed in the night in stacks of 190 frames. 

b The number of frames used in weighted shift-and-add routine described in Hor nsteinl 120061) . 
c Number of stars in initial source list. 
^ Number of stars in final source list. 

e The signal to noise ratio determined from median uncertainties of the six faintest non- variable stars detected in all the nights (SO- 14, SI -25, SO- 13, SI -68, S2-5, SI -34) with mk ~ 13.4 
mag. 



exclude source detections that occur in regions of the aver- 
age images covered by less than 50% frames that went into 
making a particular image. These regions, which are on the 
edges of the image, have relatively low signal to noise and the 
PSFs in these regions may not be well represented by the PSF 
model. We also exclude individual measurements in which 
k nown stars are blen ded with each other (i.e., sources as listed 
in lGhez et ail d2005l) as well as Sgr A* IR.) This procedure, 
in its entirety, produces 4795 detections among 131 sources, 
which range in ntx magnitude from 9.0 to 16.1 mag (see Fig- 
ure |2j. 



Photometri c zeropoints are est ablished on the basis of the 
work done by Blum et al. ( 1996a). While we share 7 stars in 
common with lBlum et all (Il996al) (IRS 16NW, 16SW, 16C, 
16NE, 29N, 29S, and 16CC), only IRS 16C (m K = 9.83 ± 
0.05 mag) is a suitable photometric referen ce source. IRS 
16SW is a known v ariable star in the K-band (lOtt et alJll999t 
iDePov et all 12004 and IRS 29N is noted as possibly vari- 
able in lHornstein et alj d2002l). IRS 16C C appears to be vari- 
able in the L-band: iBlum et alj (TT9 96a) list a re-calibrated 
value from iDepov & Sharpl dl991l) of 8.7 ±0.2 mag, while 
ISimon s & Becklin (1996) measure 10.2 ±0.2 mag. Among 




FIG. 1 . — Identification of sources from this study overlaid on a 6" x 6" region of an LGS K'-band (A D = 2. 1/xm) image from Ghez et al. 1 2006) taken on June 
30, 2005. All 131 sources are circled, but only variable sources are labeled. The location of Sgr A* is marked with an asterisk. 
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FIG. 2. — A mx histogram of the 131 stars in our sample which range in 
magnitude from 9.0 to 16.1 mag. 

the remaining sources, only IRS 16C is in the final source 
lists of all the images. Several non-variable sources (see 
§3.2.1) are used a posteriori to confirm that IRS 16C is non- 
varying. Specifically, we check for any systematic shifts in 
the zero points by examining the normalized flux densities 
(Qj = J2f Hfi^' where index i represents each star in an im- 
age of epoch j, Flux avgi is the weighted average of the flux for 
that star over all images, and N is the number of stars used) of 



the 7 least variable bright stars that are identified in all 50 im- 
ages (Sl-3, Sl-5, S2-22, S2-5, Sl-68, S0-13, and Sl-25) (see 
§3.2.1). The photometric stability ofIRS 16C is shown in Fig- 
ure [3] which plots Qj versus the observing dates. The refer- 
ence source IRS 16C appears to be stable over time, since the 
standard deviation of Qj is 0.05, which is consistent with our 
measurement uncertainty for bright stars. Increasing the num- 
ber of reference stars to 1 1 non-variable sources present in all 
frames in all 50 nights yields the same result. We therefore 
conclude that IRS 16C is non-varying to within our measure- 
ment uncertainties, and include it in our list of non-varying 
sources. Uncertainty in each of our reported relative pho- 
tometry values is initially estimated as the root mean square 
(RMS) deviation from the average of the measurements from 
the three different subset images. The RMS value is added in 
quadrature with the uncertainty in the brightness of IRS 16C 
(0.05 mag) determined from the standard deviation of the nor- 
malized flux densities Qj, As FigureHJshows, the median un- 
certainties grow from a floor of about 0.06 mag to 0.21 mag 
for the K= 16 mag sources. 
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FIG. 3. — The average normalized flux densities (Q,-, as defined in §3.1) as 
a function of time. The variations seen are consistent with the measurement 
uncertainties for bright stars and are therefore demonstrative of the stability 
of our calibrator IRS 16C. This quantity is used as a scale factor to improve 
the quality of our relative photometry. 
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FIG. 4. — The median measurement uncertainty of stars binned by magni- 
tude. The median uncertainties grow from a floor of about 0.06 mag to 0.21 
mag for the K= 16 mag sources. 

3.2. Variability 
3.2.1. Identifying Variables 

There is a wide range of methods for testing photometric 
variability and the challenge for these various approaches is to 
avoid declaring a non- vari able source variable on the basis of 
a few outlying data points dWelch & Stetson! 19 93). We there- 
fore have chosen to use the Kolmogorov-Smirnov (KS) test 
to calculate the probability that a distribution of data points 
is consistent with a model of a distribution of measurements 
for a non-variable source. This approach is less sensitive to 
outlying data points than the commonly used \ 2 test, which 
is an analysis of a single number description of how well a 
data set matches a model. In the KS test, we adopt as our 
model a non-variable light curve with gaussian-distributed un- 
certainties and we test the consistency of the measurements 
with the model. Specifically, we examine the distribution of 
X; = ux '~ ux ° r " where Flux, is the flux of a star in a im- 

age of epoch j, ctfiux] is the corresponding uncertainty, and 
Flux avg is the weighted average of the flux for that star over 
all images. The resulting KS probabilities, which have al- 
lowed values between and 1, describe how likely it is that 
a source's measurements are consistent with a non-variable 
source. Therefore variable stars, whose intensity variations 
are larger or comparable to our measurement uncertainties, 
should have very low KS probabilities. We classify a star 
as variable if it has a KS probability of less than 2.7 x 10~ 3 , 



FIG. 5. — KS test probabilities of all the stars, with low values implying 
a low probability for a star to be non-variable. We classify a star as variable 
if it has a KS probability of less than 2.7 X 10~ 3 , which is equivalent to a 
3<r cutoff value for gaussian distributed uncertainties and is marked by the 
dashed line. To ensure all our low KS probability stars are truly variable, 
we require all variable stars to have positive intrinsic variance. We mark all 
variable stars identified in this survey with hatches. 

which is the equivalent to a 3a cut for gaussian distributed 
uncertainties (see Figure 0. To ensure all our low KS prob- 
ability stars are truly variable, we also require these stars to 
have positive estimates of their intrinsic flux density variance, 

^intrinsic = Measured ~ ( 5 Y '> where me firSt term is me dispersion 

of the measured flux densities and the second term removes 
the bias introduced by the measurement uncertainties, S. At 
this point, the 7 least variable bright stars detected in all im- 
ages, used in §3.1 to define Qj, are identified. We then scale 
all our photometry by Qj in order to reduce the fluctuations 
induced by measurement errors on IRS 16C. The KS and in- 
trinsic variance tests are then repeated. Table [2] and [3] list the 
properties of the variable and non-variable stars in our sample, 
respectively, and the light curves of all variable stars and a few 
key non-variable stars are shown below in Figures [9] [10] [TT] 
[TUEJ and[]3] While there are almost certainly other variable 
stars that we have excluded, our uncertainties limit our abil- 
ity to classify more of these stars as variable, especially at the 
fainter end. 

Among the 131 stars in our sample, 15 are identified as 
photometric variables in K-band (see Figure [T]). To this we 
also add IRS 16CC, known to be variable in the L-band (see 
§3.1). Since the relative photometric uncertainties are roughly 
uniform down to a m« magnitude of ~ 14 and then grow at 
fainter magnitudes (see Figure @}, we report a frequency of 
variable stars based on the stars brighter than 14 mag. Within 
this brighter sample of 44 stars, there are 10 variable stars, 
suggesting a minimum frequency of variable stars of 23%. 
There is no evidence for radial dependence, suggesting that 
we are not limited by stellar confusion down to 14 mag. 

We compare our results to those of lOtt et all (119991) . the 
only other high spatial resolution study of the variability of 
sources in the Galactic center. Those authors give an upper 
limit of possible variable stars of approximately 50% of their 
218 sources with ntx < 13 mag over 18" x 18". While this 
variable star frequency is higher than our reported value (and 
consistent), a comparison limited to the stars in common leads 
to a number o f discrepancies. In the overlap sample of 33 
stars, lOtt et all find 2 of the stars to be variable (IRS 16SW, 
SI -3), while our sample has 6 (IRS 16SW, IRS 16NW, IRS 
29N, S2-11, S2-4, SI -21) and only IRS 16SW is in common. 
There are a number of differences between these two studies, 
including the data analysis approach used (PSF fitting vs. 
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TABLE 2 
List of Variable Stars 



Star ID 


Other ID 




Int. Var. 


P 


AR.A. 


ADec. 


Probability 


Nights 


Type 






(mag) 


(mag) 


(arcsec) 


(arcsec) 


(arcsec) 




(days) 




IRS16SW 


E23 


10.06±0.19 


0.17 


1.41 


1.04 


-0.95 


2.3E-23 


50 


Ofpe/WN9 ch 


IRS16NW 


E19 


10.09±0.09 


0.07 


1.21 


-0.01 


1.21 


3.3E-05 


49 


Ofpe/WN9 c? 


IRS29N 


E31 


10.33±0.20 


0.18 


2.15 


-1.63 


1.40 


1.1E-15 


35 


WC9 c s 


IRS16CC 


E27 


10.60±0.05 


0.01 


2.07 


1.99 


0.57 


3.0E-01 


50 


O9.5-B0.5 I c 1 


S2-11 


GEN+2.03-0.63 


11.99=1=0.13 


0.11 


2.07 


1.99 


-0.58 


4.9E-19 


49 


Late e s 


S2-4 


E28:GEN+1.46-1.49 


12.26±0.17 


0.15 


2.05 


1.45 


-1.45 


6.1E-14 


47 


BO-0.5 F« 


Sl-1 


GEN+1. 01+0.02 


13.00i0.ll 


0.08 


0.98 


0.98 


0.05 


3.6E-04 


49 


Early' 


S2-36 




13.28±0.13 


0.12 


2.08 


2.04 


0.43 


3.7E-09 


48 


Early f 


Sl-21 


E24:W7 


13.33i0.17 


0.15 


1.69 


-1.69 


0.13 


4.2E-06 


42 


09-9.5 III? CE 


Sl-12 


E2LW13 


13.82±0.18 


0.17 


1.31 


-0.85 


-1.00 


4.2E-07 


45 


OB I? c 


S2-7 


E29:GEN+1.06+1.81 


14.06±0.25 


0.21 


2.09 


0.97 


1.85 


4.2E-10 


45 


O9-B0 c 


SO-32 




14.18±0.20 


0.15 


0.81 


0.26 


0.77 


6.3E-04 


49 


Early f 


S2-58 




14.21±0.14 


0.10 


2.45 


2.17 


-1.14 


6.8E-04 


47 


Early f 


SI -45 




15.41i0.55 


0.28 


1.63 


-1.28 


1.00 


6.1E-06 


41 


Early f 


S2-65 




15.83±0.49 


0.29 


2.57 


2.37 


-1.00 


2.5E-04 


29 




SO-34 




15.85±0.40 


0.31 


0.83 


0.32 


-0.77 


4.1E-06 


26 





NOTE. — Photometry is relative to IRS 16C (m^=9.83 mag). Positions are in arcseconds offset fro m Sgr A* in 1999.56 and p is the projec ted distance. The K-S probability is equal 
to 1 for an ideal non variable sou rce, and approaches for a very variable source. O ther ID's are from Paumard et al. 1 2006) and Genzel et al. ( 2000) respectively. We classify IRS 29N 
as an early type star according to Paumard et al. 1 2006) although Fiaer et al. 1 2003) classifies it as a late type star. 

a The magnitudes are Q j corrected using 7 bright non-variable stars, and the uncertainties do not include the 5 % absolute calibration uncertainties. Comparison to other sources requires 
adding them in quadrature. Uncertainties are calculated as the standard deviation of the mean. 

Spectroscopic identification by Eisenhauer et al. 1 2005). 
c Spectroscopic identification by Paumard et al. i 2006). 
^ Spectroscopic identification by Figer et al. i 200^). 

e Spectroscopic identification by Ott 1 2003). We denote sources with clear CO or He lines as Early and Late respectively. 

Identification based on the interpretation by Genzel et al. 1 2003) of m (CO) index of Otl 1 2003) where Genzel et al. 1 2003) identify stars with m{CO) > 0.04 as late type stars and stars 
with m(CO) < 0.04 as early type stars. 
g Identified as nonvariable by Ott et al. 1 1999). 
h Identified as possibly variable by Ott et al. i 1999). 
1 IRS 16CC appears to be variable in the L-band as discussed in §3.1. 
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TABLE 3 
List of Non- variable Stars 



Star ID 


Other ID 


K 


P 


AR.A. 


ADec. 


Nights 


Type 






(mag) 


(arcsec) 


(arcsec) 


(arcsec) 


(days) 




1RS16NE 


E39 


9.00+0.05 


3.06 


2.85 


1.10 


35 


AI r\TAC Q 

Otpe/WN9 8 


IRS16C 


E20 


9.83±0.05 


1.23 


1.13 


0.50 


50 


Ofpe/WN9 cgi 
B0.5-1 I cg 


S2-17 


E34:GEN+1.27-1.87 


10.90+0.07 


2.23 


1.27 


-1.84 


35 


TD C 1 ACVU F 


C29. 1 ACF 1 


1 i nn_l_A no 


9 1 C 
Z. 1 J 


1 Q< 


111 
-i.il 


DO 


VUf^O/O^S 

WLo/y fc 


IRS29S 




1 1.31±0.06 


2.08 


-1.86 


0.93 


30 


K3 III d 


C 1 9/1 

SI -24 


Ezo.OEJN+0. /O-I.JJ 


11 A/i _i_n m 
11.64±0.0/ 


1.72 


n -71 


1 cc 


A c 

45 


(Jo-9.5 1 fc 


S2-16 


EJJ.Z9JNE1 


i i u^xn no 
ll.oD+U.Uo 


2.29 


1 m 
-1.01 


O AC 

z.05 


3J 




SI -23 


/-iT7\T n on 1 /I A 

GEN -0.90- 1.46 


1 i oi;±(i nn 
11.86+0.09 


1.73 


-0.92 


-1.46 


33 


Late eg 


S3-2 


GEN+3. 07+0.56 


12.00±0.11 


3.09 


3.03 


0.60 


31 


Early 1 
08.5-9.5 1° 


S2-6 


E30:GEN+1.60-1.36 


12.06+0.08 


2.07 


1.59 


-1.31 


50 


Sl-3 


E15:GEN+0.57+0.84 


12.10±0.06 


0.99 


0.46 


0.88 


50 


Early fh 
WN5/6 cg 


S3-5 


E40:16SE 


12.15±0.09 


3.16 


2.95 


-1.13 


21 


S2-8 


W2 


12.24±0.08 


2.16 


-1.99 


0.84 


23 


Early' 


Sl-17 


GEN+0.55-1.45 


12.51±0.08 


1.52 


0.50 


-1.44 


49 


Late eg 


Sl-4 


GEN+0.77-0.71 


12.53±0.07 


1.02 


0.77 


-0.66 


50 


Early' 


S2-19 


E36:GEN+0.53+2.27 


12.62±0.11 


2.34 


0.42 


2.30 


33 


O9-B0 I? cg 


SI -20 


GEN+0.41+1.59 


12.70±0.11 


1.66 


0.37 


1.61 


49 


Late eg 


Sl-22 


E25:W14 


12.72±().()8 


1.72 


-1.65 


-0.51 


42 


08.5-9.5 I? cg 


Sl-5 


GEN+0.43-0.96 


12.78±0.04 


0.98 


0.37 


-0.91 


50 


Late eg 


Sl-14 


E22:W10 


12.82±0.07 


1.40 


-1.37 


-0.30 


46 


08-9.5 III/I cg 


S3-6 


GEN+3. 26+0.08 


12.82±0.05 


3.22 


3.22 


0.09 


17 


Late eg 


S2-22 


GEN+2.37-0.29 


12.92+0.04 


2.33 


2.32 


-0.22 


50 


Late eg 


S2-38 




12.93±0.09 


2.12 


2.04 


0.58 


44 


Late' 


S2-31 


GEN+2.91-0.20 


13.06±0.10 


2.84 


2.83 


-0.15 


42 


Late eg 


SI -34 




1 1 in 1 A 1 A 

13.20+0.14 


1.29 


0.86 


-0.96 


50 




S2-5 


GEN+l.91-0.86 


13.32±0.05 


2.05 


1.89 


-0.80 


50 


T- 1 f 

Early 1 


SI -68 




1 J.J5±0.06 


1.97 


1.89 


-0.55 


50 




S2-21 


GEN -1.70- 1.65 


13.47±0.10 


2.36 


-1.70 


-1.65 


12 


Early 1 
Late c 


cn 1 q 
J>U-l9 


/"•fw-lA cn n /n 


1 J.4y±U.U4 


n An 

u.oy 


n <a 


n /19 
-U.4Z 


cn 


^1 9S 

k3 1 -ZJ 


flFN-i-1 6Q 66 




1 76 


1 £.5 


-u.ou 


co 




S2-26 






2.56 


0.69 


2.47 


20 


Late e 


SO- 15 


E16:W5 


13.70±0.07 


0.98 


-0.93 


0.29 


48 


09-9.5 V c 


SO-14 


E14:W9 


13.72±0.08 


0.83 


-0.78 


-0.27 


50 


09.5-B2 V c 


Sl-19 


GEN+0.38-1.58 


13.82±0.12 


1.62 


0.36 


-1.58 


46 


Early' 


S2-2 


GEN-0.54+2.00 


14.07±0.11 


2.12 


-0.59 


2.03 


41 


Late' 


SO-2 


E1:S2 


14.16±0.08 


0.12 


-0.07 


0.10 


33 


BO-2 V b 


Sl-8 


E18:W11 


14.19±0.11 


1.08 


-0.67 


-0.85 


49 


OB c 


Sl-15 


W4 


14.21±0.10 


1.46 


-1.37 


0.52 


47 


Late' 


SO-6 


S10 


14.26±0.09 


0.39 


0.07 


-0.38 


49 


Late' 


Sl-49 




14.26+0.13 


1.66 


-1.65 


0.15 


23 




Sl-13 


W12 


14.27±0.12 


1.42 


-1.10 


-0.90 


46 


Early' 


S2-47 




14.29±0.08 


2.26 


2.20 


-0.49 


48 


Early' 


SO-9 


Sll 


14.31±0.08 


0.55 


0.14 


-0.53 


49 


Early' 


SO- 12 


W6 


14.38±0.06 


0.68 


-0.57 


0.37 


49 


Late' 


S2-3 


W15 


14 48+0 09 


2.09 


-1.54 


-1.41 


23 


Late' 


c.a a 


fi a-o.o 


1 A AQ-\-C\ 1 1 


a ii 


n 39 


1 Q 
-U. ly 


Act 


RA 9 vb 

DU-Z V 




FA-C/1 


14.jU_tU. 14 


n 9^ 
U.Zj 


n 99 
U.Zz 


n 1 1 

0. 19 


1 Q 


DA 9 \/b 
DU-Z V 


c.9 7C, 




1 A S94-0 1 9 


9 78 
Z. / 


9 £.5 

Z.OJ 


85 


40 




Sz-69 




14.j/±U.i4 


9 A/1 
Z.t)4 


n ni 
-0.91 


9 A O 
Z.45 


1 1 

1 J 


Early' 


S3-4 




14.61+0.20 


3.14 


3.10 


-0.47 


20 


Early 1 


S0-1 


E4:S1 


14.67+0.11 


0.14 


-0.11 


-0.09 


49 


BO-2 V 


S2-23 




I4.72±0.l3 


2.43 


1.64 


1.80 


39 


Late' 


Sl-55 




14. 5U±U. 39 


i An 
1.69 


1 CO 
1.90 


n en 
0.59 


A 1 

41 




c 1 cn 




1 A QTJ_n QT 

14.5Z±U.-? / 


1 A7 
1.0/ 


1 CI 
1.91 


n 79 
0. /z 


A 1 

41 




Sl-52 






1.66 


-0.02 


1.66 


42 




Sl-10 


W8 


14 so 4-0 12 


1.15 


-1.15 


-0.04 


42 




9 


Fi7-nFw n A6 1 ai 

IL1 / .UIJ.1N-U.UO- 1.U1 




1 OA 


-U.U9 


1 no 
-1 .uu 


4S 


i f 
Early 


Sl-33 






1.25 


-1.24 


-0.07 


40 


Sl-58 




1 c 04+0 35 


1.77 


-1.48 


0.98 


37 




Sl-51 




1 5 0S4-0 1 'S 


1.66 


-1.65 


-0.20 


39 


F 1 ' 


S2-86 




1 S 0S4-0 73 


2.99 


2.68 


-1.33 


24 








1 J .UuxU. jj 


1 1 9 


3 Afi 


69 
-u.uz 


1 Q 


if 
Early 


S2-30 




15.13+0.19 


2.88 


2.88 


0.00 


23 


Sl-18 




15.14+0.15 


1.66 


-0.73 


1.49 


39 


Early' 


S0-5 


E9:S9 


15.17+0.16 


0.36 


0.18 


-0.31 


43 


BO-2 V b 


S0-31 


E13 


15.20+0.23 


0.66 


0.49 


0.45 


40 


B V c 


S2-34 




15.21+0.22 


2.04 


1.79 


0.98 


43 




SO-26 


E8:S5 


15.27+0.22 


0.39 


0.36 


0.16 


40 


B4-9 V b 


SI -44 




15.28+0.41 


1.61 


0.26 


1.59 


40 




S3-16 




15.30+0.20 


3.15 


3.02 


-0.88 


14 


Late' 


S2-82 




15.30+0.29 


2.88 


2.87 


0.08 


27 


Late' 
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List of Non- variable Stars 



S2-12 




15.33±0.17 


2.07 


1.68 


1.22 


38 


Late f 


Sl-32 




15.33±0.13 


1.13 


-0.93 


-0.64 


42 




Sl-39 




15.35±0.16 


1.45 


-0.53 


-1.35 


30 


Early f 


SO- 11 


E12:S7 


15.36±0.13 


0.53 


0.53 


-0.01 


40 


B V c 


SO- 18 


S18 


15.36±0.13 


0.43 


-0.09 


-0.42 


45 




Sl-35 




15.36±0.19 


1.27 


-1.24 


-0.25 


38 




S2-63 




15.39±0.28 


2.56 


-0.69 


2.47 


16 


Early f 


SI -54 




15.41±0.23 


1.68 


-1.53 


0.70 


36 




SI -62 




15.41±0.31 


1.82 


0.51 


1.74 


37 




Sl-53 




15.44±0.15 


1.68 


1.68 


-0.09 


30 




S2-61 




15.46±0.17 


2.54 


2.46 


-0.64 


38 




S2-46 




15.47±0.30 


2.18 


2.08 


-0.64 


38 




S2-73 




15.50±0.32 


2.72 


2.21 


-1.58 


29 




Sl-6 




15.54±0.13 


1.15 


-0.91 


0.71 


34 


Early* 


SI -48 




15.54±0.23 


1.62 


-0.60 


-1.50 


22 




SO-7 


E11:S6 


15.55±0.27 


0.45 


0.44 


0.11 


35 


B V c 


SO- 19 


E5 


15.56±0.19 


0.19 


-0.09 


0.17 


28 


B4-9 V b 


SI -64 




15.57±0.36 


1.91 


0.60 


1.81 


34 




S2-80 




15.57±0.27 


2.86 


2.22 


1.80 


20 


Early* 


S2-40 




15.58±0.18 


2.15 


1.68 


1.34 


35 


Early* 


S2-42 




15.61±0.28 


2.11 


0.41 


2.07 


26 




SI -27 




15.62±0.27 


1.09 


-1.07 


0.24 


39 




SI -26 




15.62±0.13 


1.02 


-0.95 


0.38 


40 




Sl-37 




15.63±0.16 


1.42 


-1.34 


0.47 


36 




SI -47 




15.67±0.23 


1.63 


-1.57 


0.45 


34 




SO- 16 


E2 


15.68±0.19 


0.08 


0.05 


0.06 


18 


B4-9 V b 


Sl-59 




15.70±0.25 


1.86 


0.01 


1.86 


27 




Sl-31 


GEN-0.91+0.44 


15.70±0.18 


1.14 


-0.99 


0.57 


31 




SO-29 




15.72±0.50 


0.54 


0.25 


-0.48 


21 




S2-83 




15.74±0.19 


2.94 


2.87 


-0.63 


16 




SO-27 




15.74±0.18 


0.54 


0.13 


0.52 


32 




Sl-36 




15.75±0.27 


1.36 


-0.67 


-1.18 


30 




S2-37 




15.77±0.31 


2.11 


0.04 


2.11 


29 




SO-8 


E7 


15.79±0.14 


0.47 


-0.32 


0.34 


29 


B4-9 V b 


Sl-7 




15.81±0.14 


1.12 


-1.00 


-0.50 


34 




SI -65 




15.82±0.15 


1.93 


1.43 


1.29 


31 


Early* 


SO-28 


S19 


15.85±0.18 


0.61 


-0.18 


-0.58 


30 




S2-64 




15.85±0.59 


2.56 


2.54 


0.31 


16 




SO-36 




15.85±0.28 


1.03 


-0.60 


-0.84 


29 




SO-20 


E3 


15.86±0.20 


0.21 


-0.18 


-0.10 


33 


B4-9 V b 


SI -40 




15.95±0.23 


1.50 


-1.36 


-0.64 


26 




Sl-61 




15.99±0.38 


1.76 


-1.43 


-1.03 


21 




S2-52 




16.02±0.20 


2.37 


2.37 


-0.07 


21 




SI -42 




16.13±0.22 


1.60 


0.94 


1.29 


19 





NOTE. — See notes from table^ 
3 The magnitudes are Q } corrected using 7 bright non- variable stars, and the uncer- 
tainties do not include the 5 % calibration uncertainties. Comparison to other sources 
requires adding them in quadrature. 

b_h Type from fPaumard et alj O00^ . [Hge7eTaT] f2003h . lEisenhauer eTaT] f2Q0l . and 
[OtH 420031) as described in table[T] 
g Identified as nonvariable by Ott et al. 1 1999). 
k Identified as possibly variable by lOtt et^T]|199^ . 

1 This is our main calibration star and is included in this table only for completeness. 

The 7 least variable bright stars detected in all images used for scaling of the photom- 
etry in order to reduce the fluctuations induced by measurement errors on IRS 16C (see 
§3.1). 



Stellar Variability in Galactic Center 



9 



15 




I , , i I i i i I , , , I , , , I , , , I i i i I i 

20 40 60 80 100 120 

period (days) 

FIG. 6. — D etermination of IR S 16SW's period using Dworetsky's string 
length method Dworetsky 1 1983). Phi is the string length where shorter 
lengths and therefore smaller values of Phi indicate a periodic signal. Note 
that the other harmonics are also found at multiples of the original period. 

aperture photometry), the time baseline (10 years vs. 5 years), 
and the angular resolution (0"05 vs. 0"13 ). Since our study 
covers twice the time baseline, we can pick out variations on 
longer time scales, which helps explain our additional vari- 
ables. Also, the high stellar crowding makes the area we 
observe the most uncertain region for the lower resolution 
lOtt et alj study. Only one star, Sl-3, in our non-variable sam- 
ple is identified as variable bv lOtt et ail and it is very close to 
their threshold for variability. 

3.2.2. Variability Characterization and Periodicity Search 

We attempt to characterize the minimum time-scale for 
variation by searching for daily and monthly variability using 
KS tests similar to the test for variability in §3.2.1 where we 
adopt as our model a non-variable light curve with gaussian- 
distributed uncertainties. For the daily variations, we group 
the consecutive nights in pairs and examine the distribution 
of the pairs in sets i for each star of X,- = F '" x ' F/ " A * where 

FluXj and Fluxf, are the fluxes of stars in images of consec- 
utive epochs j and k respectively, UFiuxj and <7f/„ A( are the 
corresponding uncertainties. For the monthly variation, all 
measurements made within days of each other are averaged 
together and the pairs separated by one month are examined 
with the same KS test. The only stars showing daily vari- 
ability in excess of 3<r is IRS 16SW, and the stars showing 
monthly variations are IRS 16SW and S2-36. 

The light curves of the variable stars are searched for peri- 
odicities using three different m ethods. First, the Lomb and 
Scargle periodo gram technique (lLomblll976t IScargldll982t 
iPress et all 19921) . which fits fourier components to the data 
points, is applied and is expected to yield a larger power 
spectral density at intrinsic harmonics of a data set in which 
there is a periodic signal. Second, Dworetsky's (1983 ) string 
length method, a variant of the Lafler-Kinman method ( 1965b . 
phases the data for every possible period and then sums 
over the total separation between points in phase space, with 
the best period and its alia ses co rresponding to the smallest 
lengths. Third, Stetson's d 19961) string length technique is 
similar to Dworetsky's, but also weights these lengths by their 
uncertainties and how close in phase the points are. Our crite- 
rion for considering a star periodic is that it show similar pe- 
riods from all three techniques. The periodicity search shows 
only one periodic star: IRS 16SW. We find a photome tric pe- 
riod of 9.724 ± 0.001 days , which is consistent with OttetaL 
( 1999), De Pov et al.1 d2004l) and the reanalysis of the btt et al 



data in iMartins et al.l d2006l) . Figure [6] shows the determina- 
tion of IRS 16SW's period using Dworetsky's (1983) string 
length algorithm, showing the 9.724 day period and its other 
harmonics at multiples of its period. The top panel in Figure|7] 
shows the phased light curve of IRS 16SW at 9.724 days and 
depicts a clearly periodic signal with an amplitude of ~ 0.55 
mag. The phased light curve of IRS 16SW is asymmetric (see 
Figure[7]l with a rise-time that is ~ 1.6 times longer than the 
fall-time. 

4. DISCUSSION 

The 16 variable stars identified in this study cover a wide 
variety of different types of stars, as we only limited our 
search by location and brightness. As Figure [8] shows, based 
on the K magnitudes alone, this sample is expected to con- 
tain early-type (O & B) main sequence stars, late-type (K & 
M) giant stars, and most types of supergiants. Fortunately, 
all but two of the variable stars have spectral classifications 
(see column 10 in Table|2]i. While 9 of the va riable stars are 
secur e ly identified from spec t roscopic work (iPaumard et al] 
I2006t lEisenhauer etail 120051; iFiger et all 120031; lOttl 1200*31) 7 
an additional 5 stars are classified on the bas i s of narrow- 
band photometry of CO absorption (Ott 20031 iGenzel et all 
120031) . In summary, four stars are LBV candidates (IRS 
16SW, 16NW, 16NE, & 16C; see §4.1) one is a WC9 (IRS 
29N; see §4.2), four are OB supergiants (S2-4, Sl-12, S2-7, 
& IRS 16CC; see §4.3), one is an O giant (SI -21), five more 
are classified as some sort of early-type star from narrowband 
filter photometric measurements (Sl-1, S2-36, SO-32, S2-58, 
SI -45), and one is classified as a late-type star from narrow- 
band filter photometric measurements (S2-1 1; see §4.3). The 
variability of stars in each spectral classification is discussed 
in turn below, along with a discussion of the possibility of 
external agents causing variability. 

4.1. Ofpe/WN9 Stars 
4.1.1. Luminous Blue Variable Candidates 

Four stars in our sample are Ofpe/WN9 stars and have been 
previously classified as candidate LBVs (IRS 16NE, 16C, 
16SW, and 16NW) based on their bright luminosity, their nar- 
row e mission lines, and t heir proximity and sim i larity to IRS 
34W (IClarketalJl200l IPaumard et all l2004bt rTrippe et alJ 
120061) 5 . In our observations, both IRS 16NE and 16C are 
non-variable over a ten-year time frame to within our un- 
certainties. The other two LBV candidates (IRS 16NW and 
16SW) show variability, but not the characteristic LBV erup- 
tions, which in the context of this study are the AM,, ~ 1 - 2 
mag events oc curring every 10-40 years and l asting as long as 
several years (Humphrevs & Davidson 1994). IRS 16SW has 
periodic variability that is explained by an eclipsing binary 
system (see §4.1.2 below) and IRS 16NW has an overall flat 
light curve with a decrease in brightness (Amx — 0.2) mag 
between 1997 and 1999 (see Figure [9]). The apparent dim- 
ming of IRS 16NW can be explained by ejected circumstellar 
material obscuring the star, with an amplitude that is smaller 
than is characteristic of a typical LBV. These variations do 
not require it to be an LBV, just that it has strong stochastic 

5 The classification of IRS 34W as an LBV is based on its bright lu- 
minosity, narrow emis s ion lines, along wi th a multi-year obscuration event 
iPaumard et al. 2004b; Trippe et al. 2006). However, more recent studies 
have cast doubt on IRS 34W's categorization as an LBV since it lacks spectro- 
scopic variability and si nce the eruption ev ent responsible for the obscuration 
event was not observedl Trippe et al. 2006). 
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mass stars in an eclipsing binary star system. 
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winds. While none of these stars shows the classic character- 
istics of LB Vs, our time baseline is too short to rule them out 
as LBVs as they may be in a quiescent phase. Nonetheless, 
our observations do not demand the complication of multiple 
recent star formation events that would be suggested with the 
presence of both LBVs and WR stars. 

4.1.2. Asymmetric Periodic Light Variations in IRS 16SW: Tidal 
Deformation ? 

The periodic variation in IRS 16SW has recently been 
attributed to either an equal mass contact ec li psing-binary 
or a massive pulsating s tar dOtt et al.l 119991; iDePov et al.l 
12004 iMartins et ail 120061) . although t he measurement of a 
spectroscopic radial velocity period by Mar tins et al.l (12006b 
strongly suggests it is an eclipsing binary star system. Also, 
a re-analysis of the data suggesting a pulsating star now 
agrees with IRS 16 SW probably being an eclipsing binary 
(Peeples et al. 2006). The asymmetries that we observe in the 
phased light curve of IRS 16SW are difficult to explain in the 
context of an eclipsing binary star system. Figure|7]shows the 
properly phased light curve (19.448 ±0.002 days), in which 
asymmetries in its rise and fall-times are still readily detected. 
The asymmetry is remarkably similar for both halves of the 
phased light curve. While this was not detected in earlier pho- 
tometric studies 6 , our study is likely more sensitive to small 
photometric variations due to our higher angular resolution 
(see §3.2.1). The observed asymmetry has been sustained 
over 10 years; if the asymmetry has been changing over time, 
it would show up as a dispersion in the vertical placements of 
the points in the phase diagram that is much larger than what 
we observe. In figure [7] we explicitly show that the asymme- 
tries are the same during the first half and second half of the 
data-set and we see no period drifts, with the period of the 
first half and second half not differing at the 3er confidence 
level. Magnetic hot spots can explain ligh t curve asymmetries 
(Diurasevic et alj2000tlCohen et alj|2004l) . however, they can 
not maintain this asymmetry over long periods of time. Fur- 
thermore, the similarity in the asymmetry between the first 
and second half would require the spots to be the same on 
both stars. Likewise, heating in the contact region of the bi- 
nary or any third light in this region produces a light curve 
that is a mirror re flection between the 1st and 2nd half (see 
Moffat et al. 2004) and is therefore inconsistent with our ob- 
servations. Other heating mechanisms such as irradiation ef- 
fects are unlikely to be the c ause of the a symmetry as they 
barely change the light curve (Bauer 2005). Given that this is 
suspected to be a contact binary, we suggest that the asymmet- 
ric light curve may be due to tidal deformations caused by the 
proximity of the stars in asynchronous orbits. Since these two 
stars are equal in mass ( ~ 50M(^), equal in r adius (66Rq) and 
appear to be in contact (IMartins et al. 2006), they fall within 
the tidal radius of ~ 163/?©, and are likely tidally deformed. 
In order to produce the asymmetric light curve, the rotation 
of the stars needs to be asynchronous with their orbital pe- 
riods so that the rotational inertia of the stars prevents their 
tidal bulges from being aligned with the line joining the stars' 
centers of mass. Two stars in such close contact will eventu- 
ally synchronize their orbital period with their rotation period. 
The synchronization time scale of stars with convective cores 
and radiative envelopes is longer than for stars with convec- 



6 Our light curve is similar to the light curve presented by Ott et al. 1 1999), 
althoug h the asymmetry wa s not explicitly reported and the reanalysis of the 
data by i Martins et al. 2006) does not show the asymmetry. 



tive envelopes, albeit more difficult to calculate (Zahn] ll977l) . 
We therefore calculate the convective envelope synchroniza- 
tion time as a minimum time for s ynchronization based on 
formalism developed by Zahn (1 1977b and find a minimum age 
of ~ 6 x 10 s yrs. This time scale is far larger than the lifetime 
of stars as massive as IRS 16SW and asynchronous orbits are 
therefore acceptable. The equality of the two halves of the 
light curve implies that the rotation rate of the two stars is 
very similar. Since these two stars are identical in every other 
respect, this equality is not surprising. It is therefore possible 
that the asymmetric light curve is due to tidal deformations. 

4. 1 .3. Eclipsing Binary Fraction 

Stars close to the Galactic center are expected to have a 
higher fraction of ellipsoidal and eclipsing variable binaries 
than the stars in the solar neighborhood. However, IRS 16SW 
is the only eclipsing binary star detected in our sample. The 
stars' orbits may be smaller due to hardening by encounters 
with other stars producing tightly bound binaries, such that 
eclipses would be more likely. In addition, collisions and 
tidal capture produce binar ies, making ellips oidal variations 
or eclipses more probable (iTamura et alj|1999l) . Of the 164 
Galactic O s t ars in clusters or associations in the sample by 
iMason et al.1 (fl998). 50 are confirmed as spectroscopic bina- 
ries, 40 are unconfirmed spectroscopic binaries, 4 are con- 
firmed eclipsing variables, and 14 are either ellipsoidal or 
eclipsing binaries. This yields local rates of eclipsing O star 
binaries between 2% and 11%. Our sample contains 1 1 spec- 
troscopically confirmed O stars, 4 Ofpe/WN9 stars, and pos- 
sibly more unconfirmed. If the Galactic center fraction of 
eclipsing binaries is similar to the cluster results, we would 
expect on the order of one eclipsing binary. Our detection 
of one eclipsing variable star suggests that the frequency of 
eclipsing binaries is not significantly increased at the Galactic 
center over the local neighborhood. 

4.2. Late-type WC Wolf-Rayet Stars: Variations Associated 
with a Wind Colliding Binary 

Three stars in our sample are spectroscopica lly identified as 
WC s tars (IRS 29N, IRS 16SW-E and S2-16) dPaumard et alj 
2006), and all three a re dust producers as evidenced by thei r 
red colors (K-L ~ 3; iBlum et al.|[~1996al IWright et al.l 2006). 
In this study, only IRS 29N is variable (see Fig[l0]>. IRS 29N's 
intensity shows a gradual drop and then rise in brightness of 
Am.K — 0.7 mag over a time scale of « 5 years. Its light 
curve is similar to the variations s een in WC stars elsewhere i n 
the Galaxy (e.g. compilation by Ivan der Hucht et al.l l2001b); 
these sources are thought to be variable due to periodic or 
episodic dust production in the wind collision zone of long 
period eccentric binary star systems (lOOOd < P < lOOOOd) 
during periastron passage. When the dust forms, the star ex- 
hibits a rising infrared flux followed by fading emission when 
dust fo rmation stops and dust grains are dispersed by stellar 
wind s dMoffat et al.lll987t IWhite & Beckerill995t I Veen et al.l 
[l998l;IWilliams & van der Huchtl 19921 l200ol) . Currently, only 
seven WC stars have been observed to produce dust episodi- 
cally, all of which are c onfirmed or suspected massive binaries 
with elliptical orbits (|van der Huchtl 12001 at IWilliams et al.1 
2005] iLefevre et al.l 120051) . Two of these are known to ex- 
hibit p inwheel nebulae, a tell-tale sign of win d-colliding bina- 
ries dTuthill et al.ll999tlMonnier et al.ll999l) . The time-scales 
and magnitude of the photometric variability of IRS 29N is 
consistent with it being a wind-colliding binary; we therefore 
conclude that IRS 29N is likely to be a wind-colliding binary. 
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FIG. 9. — Light curves of the four stars in our sample that are classified as candidate LBVs (IRS 16SW, 16NW, 16NE, 16C). The top two panels denote the two 
variable stars, while the bottom two panels are non-variable stars. The given intrinsic variance is in magnitudes. 
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FIG. 10. — Light curves of Wolf-Rayet stars of type WC in our sample. The first panel is the variable star IRS 29N, which is probably a wind-colliding binary, 
while the other two are non- variable stars. The given intrinsic variance is expressed in magnitudes. 
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FIG. 1 1 . — Light curves of variable OB stars in our sample. IRS 16CC is not variable in K band, but shows variability at L (see §3.1). We propose that S2-4 
and S 1-12 are variable due to ejection of circumstellar material obscuring the stars, while S2-7 is probably a Be star with the formation of an equatorial disk 
obscuring the star. The given intrinsic variance is expressed in magnitudes. 



14 



Rafelski, Ghez, Hornstein, Lu, & Morris 



■ Variable 

■ SpT: Early 



1996 1998 2000 2002 2004 2006 

Date (year) 



1996 1998 2000 2002 2004 2006 

Dote (year) 



4.0 -<> § 



1998 2000 2002 

Date (year) 



2004 2006 



1998 2000 2002 2004 2006 

Date (year) 



14.5 -Variable 

- SpT; Early 



<> ii 



<> <> 



S 15.5 -"" 
E -<> 



U> <b 

<> ^ 



1996 1998 2000 2002 2004 2006 

Date (yeor) 



FIG. 12. — Light curves of early-type variable stars in our sample. S2-36 is variable on a monthly basis and its variation is therefore probably intrinsic. The 
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4.3. Comments on Other Stars 

Our sample includes at least 10 other variable young mas- 
sive stars in addition to those discussed in §4.1 & §4.2. Five 
are spectroscopically identified O B stars (IRS 16CC 7 , S2-4, 

51- 21, Sl-12, S2-7; see FigureUD dPaumard et al.ll2006h . and 
five are early-type stars classified as on the basis of narrow- 
band filter photometric meas urements (Sl-1, S2 - 36, SO-32 , 

52- 58, Sl-45; see Figure E) dGenzel et al.ll200l IOttll200l . 
The majority of these variables are likely associated with mass 
loss, although interstellar extinction could also play a role, as 
discussed below. In particular, the two K-band variable OB 
supergiants each show a dip (0.3 - 0.9 mag) in their bright- 
ness that lasts for 1-6 years. OB supergiants have stochas- 
tic winds with high mass loss rates (0.2-20 x lO~ 6 M /yr) 
(Massey 2003); it is therefore likely the large variations seen 
are due to ejection of circumstellar material. Another exam- 
ple of an OB star with variations potentially due to mass loss 
is S2-7, which is either luminosity class III or V based on 
its assumed distance. S2-7 shows a decrease in luminosity 
between 2000 and 2005 with Am K ~ 0.5-1.0 mag. This is 
reminiscent of Be stars, which sometimes show fading events 
due to the for mation of an equatorial disk on time scales of 
several years (lMennickentetal.1119941: iPavlovski et al.l [19971: 
iPercv & Bakos 2001). The remaining young variable stars in 
our sample have variations that are more difficult to charac- 
terize, although S2-36 seems to have significant variations on 
monthly time-scales (see §3.2.2). Nonetheless, it is likely that 
mass loss also plays a central role in generating the observed 
variations. 

4.4. Interstellar Material 

4.4.1. Apparent Variations Caused by Stellar Motion Through the 
Line of Sight Extinction 

Periods of reduced luminosity in stars can also be due to 
external effects such as obscuration by foreground interstel- 
lar matter. The central parsec has several gas patches that 
are a few arcseconds wide and that can cause l ocal extinction 
enhanc ements of ~ 1 magnitude in K-band dPaumard et al.1 
2004a). As our FOV is only 5"x5", these patches would 
cover a sizable region and many neighboring stars would 
show similar variations. Interstellar material closer to Earth 
in the line of sight is unlikely, as i ts proximity woul d make 
the clouds even larger in projection (Trip pe et al.l2 006). Since 
neighboring stars do not experience similar effects of obscura- 
tion, the large scale structure in the interstellar medium (ISM) 
is unlikely to be the cause. 

In some cases it is possible that the variable obscuration is 
caused by the relative motion of the foreground high density 
streamers and the background stars, such as are observed in 
the L-band associated with the Northern A rm dClenet et al.l 
12004 iGhez et al.ll2005bt iMuzic et al.ll2006l) . In this area we 
see unresolved streamers that are small in projected width 
(< 80 mas). These streamers may be due to shocks heating 
the neighboring dust with the streamers tracin g thin shells 
of co mpressed gas from one or several shocks dClenet et al.l 
2004). These streamers may cause dips in our light curves 
due to small-scale structure obscuration in the line of sight. 
This causes photometric variability due to the relative lateral 
motion between the absorbing feature and the star. The pro- 

7 IRS 16CC is not identified as variable in this survey, although re- 
ported differe nces in L-band magnitudes from previous studies suggest it is 
variable IDepov & Sharrjfl99lL Simons & Becklin 1996; Blu m et al.|[T996al : 
Wright et al. 2006) (see §3.1). 



jected implied width of the small-scale structure is approxi- 
mately ~ 10-40 mas assuming a projected stellar velocity of 
~ 5 mas/yr which is typical of stars at projected distances of 
~ 2" from Sgr A*. The three stars whose variability is most 
likely ascribable to these thin high density streamers are S2- 
11, Sl -45, S2-58 . The most clear case is the late-type star 
S2-11 (Ott 2003), which shows an interval of reduced lumi- 
nosity between 2001 and 2005 with Am^ ~ 0.3 mag and is 
otherwise constant over our time frame (se e FigureQj) . Using 
a Galactic center distanc e of r = 8.0 kpc (lReidll 19931) and an 
extinction Ak=3. 3 mag dBlum et al.ll 1 996al) . we determine its 
spectral type based on luminosity and late-type classification 
as M3-5 III. Stars with spectral type M5 and luminosity class 
III are generally classified as asymptotic giant branch stars 
(AGB), but the variations observed are not t ypical of AGB 
stars which have periods between 0.5 - 1.5 yrs (Habing 1996). 
This s tar is likely obscu red by dust given its red color (K-L 
~ 2; IWright et al]|2006l) . It is located in the middle of the 
Northern Arm (see Figure fT4l > and its reduced luminosity is 
likely due to a high density streamer in the line of sight. The 
two other stars whose variability can probably be attributed 
to high density streamers have long term variations over ten 
years; S 1 -45 appears to brighten by Am^ — 0.5 - 1 .0 mag and 
S2-58 appears to dim by Amj- ~ 0.2 mag (see Figure [12). It 
is also possible that the dips in the light curves of stars specu- 
lated to be variable due to high stellar winds such as S2-4, Sl- 
12, S2-7, and S2-36 are actually variable due to obscuration 
in the line of sight. Measurements at multiple wavelengths 
throughout future variations would help to establish the role 
of variable extinction in the observed K-band variations. Fur- 
thermore, foreground material would be polarized due to the 
magnetic fields at the Galactic center, and therefore polariza- 
tion variations of stars would provide a test of the hypothesis 
that the relative motion of streamers and stars are responsible 
for stellar intensity variations. 

4.4.2. Variability of Stars Near Closest Approach 

We detect no variability in the 7 central arcsecond sources 
that have known 3-dimensional orbits (S0-1, SO-2, SO-4, S0- 
5, SO-16, SO-19, SO-20) (see Figure [B). The three fainter 
stars (SO-16, SO-19, SO-20) have missing measurements that 
are due to insufficient image sensitivity, although they are de- 
tected in highe r signal-to-noise i mages made from multiple 
nights of data (Ghe z et al.ll2005l) . The photometry of these 
stars constrains the properties of a cold, geometrically-thin 
inactive accretion disk around Sgr A*, since in the presence 
of such a disk we would expect to see nearby stars signif- 
icantly flaring in the NIR as they passed through and inter- 
acted with the disk, and eclipsed at other times. When a 
nearby star approaches such a disk we would see enhanced 
NIR flux from reprocessed UV and optical starlight incident 
on the disk (which we call a flare). Also, we would expect 
the disk to eclipse the star, reducing the flux from the star in 
varying amounts depending on the properties of the disk. The 
time-scales vary based on the geometry of the disk but are on 
the order o f a year and months for the e c lipses and flares, re - 
spectively dNavakshin & Sunvaevll200l ICuadra et al]|2003l) . 
An optically thin disk may not fully eclipse stars, and gaps in 
our observations would allow different geometries of the disk 
to account for any one st ar not showing ecl ipses or flares as 
was calculated for SO-2 (Cua dra et al.ll200 3). However, with 
the ensemble of stars that we have monitored, the effects of a 
disk with any orientation should be evident. We constrain the 
NIR optical depth to be < 0. 1 . If we assume the NIR standard 
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FIG. 13. — Light curves of a late-type star, and two unknown type stars in our sample. The late-type star, S2-1 1, is probably an AGB star of spectral type M3-5 
III, while we are unable to classify the two stars of unknown spectral type. The given intrinsic variance is expressed in magnitudes. 
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FIG. 14. — Identification of the source S2-11 from this study overlaid on 
a 5" x 5" region of an LGS L'-band (A„ = 3.8/mi) image from G hez et aT] 
1 2006) taken on June 30, 2005. The star S2-1 1 is situated in the middle of the 
Northern Arm and it is therefore likely that the interval of reduced luminosity 
is due to the stellar motion through a high density streamer in the line of sight. 
The location of Sgr A* is marked with an asterisk. 

interstellar dust opacity at 2.2fim where the dust extinction is 
ap proximately ~5x 10~ 22 cm 2 per hydrogen atom (see Fig. 2 
in Voshchinnikov et al. 2003), we find the column density of 
the disc to be ~ 2 x 10 cm~ 2 . However, dust grains may be 
larger in size or non-e xistant in the disc (ICuadra et al.ll2Q03k 
iNavakshin et al.ll2004l) . and such constraints should therefore 
be approached cautiously. Regardless, the lack of observed 



flares or eclipses in the 7 central arcsecond sources that have 
known 3 -dimensional orbits puts such severe constrains on the 
density and size of any possible disc disk around Sgr A* that 
such a disk is unlikely to exist. 

5. SUMMARY 

We use ten years of diffraction-limited K-band speckle data 
to determine the photometric stellar variability in the central 
5" x 5" of our Galaxy. Within this study's limiting mag- 
nitude of mk < 16 mag, we find 15 K-band variable stars 
out of 131 well-sampled stars. Among 46 stars brighter than 
niK < 14 mag with uniform photometric uncertainties, there 
are 10 variable stars, suggesting a minimum variable star fre- 
quency of 23%. We find one periodic star, IRS 16S W, with a 
period of P=l 9.448 ± 0.002 days, in agreement with lOtt et al.1 
(119991) . iDePov et all (|2004|) . and iMartins et all (120061) . Our 
data are consistent with an eclipsing binary and show a rise- 
time that is ~ 1 .6 times longer than the fall time and we sug- 
gest that the asymmetric light curve results from tidal defor- 
mations of the two stars in the presence of asynchronous ro- 
tation. We expect to see on the order of one eclipsing bi- 
nary in our sample for conditions similar to the rest of the 
Galaxy, suggesting that the frequency of eclipsing binaries is 
not significantly increased at the Galactic center over the local 
neighborhood. We identify IRS 29N as a wind colliding bi- 
nary based on its light curve and spectral classification. This 
rare object warrants further investigation to confirm its binary 
nature. None of the IRS 16 stars shows the classic eruptive 
events of LBVs, although our time baseline is too short to 
rule them out as LBVs. Among the remaining variable early- 
type stars in our sample, 3 exhibit large variations on time- 
scales of a year, which are either due to obscuration from mass 
loss events or from line of sight extinction variations. Three 
more stars in our sample exhibit long term variations of ~5- 
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FIG. 1 5 . — Light curves of 7 central arcsecond sources that have known 3-dimensional orbits. The X's mark observations that are rejected due to confusion with 
other sources, and the boxes are areas where the nightly images have missing measurements due to insufficient image sensitivity and correspond to detections in 
monthly averaged images. None of these stars are variable, which is evidence against the presence of a dusty disk around Sgr A* . The given intrinsic variance is 
expressed in magnitudes. 



10 yrs probably due to line of sight extinction variations due 
to high density streamers. Seven stars in the central arcsecond 
do not show photometric variations indicative of a cold geo- 
metrically thin inactive accretion disk which puts such severe 
constraints on the density and size of any possible disk around 
Sgr A* that such a disk is unlikely to exist. 
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